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Abstract 

New results, namely the independent determination of the deuterium abundance 
in several quasar absorption systems, and the complementary determination of the 
cosmological baryon density by observations of anisotropics in the cosmic microwave 
background (CMB), allow for a reevaluation of the constraints on the relativistic 
particle content of the universe at primordial nucleosynthesis. Expressed in terms 
of the neutrino energy density, we find 1.7 < N u < 3.5 (95% CL). In particu- 
lar, we show that phenomenological four neutrino models including a sterile state 
(not participating in SU(2)l x U(1)y interactions) unavoidably thermalize a fourth 
neutrino, and are highly disfavored in the standard minimal model of primordial 
nucleosynthesis, if the systematic uncertainty in the primordial helium abundance 
is small. We describe plausible extensions of the minimal model which evade this 
constraint. 



1 Introduction 



With the simplifying conditions of isotropy, homogeneity, thermal equilib- 
rium, and the particle content of the standard model of particle physics, big 
bang nucleosynthesis (BBN) is a one parameter model dependent only on the 
baryon-to-photon ratio, rj, at that epoch. The appeal of this simple yet suc- 
cessful standard model [1] has motivated a predictive ability, for example, in 
the number of leptons in the standard model of particle physics [2] . 

The light nuclides D, 3 He, 4 He and 7 Li are produced in measurable quantities 
in the first three minutes of the standard cosmology. Considerable attention 
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has been devoted to the analysis of uncertainties in the predicted abundances 
of the light elements and their consistency with the observed light element 
abundances [3,4,5,6]. Though systematic uncertainties likely dominate the he- 
lium abundance measurements [6,7,8] and may be present in observations of 
the deuterium-to- hydrogen ratio D/H in high-redshift quasar systems [9], and 
7 Li may be partially depleted by stellar processes [10], standard BBN is re- 
markably successful predictor for the abundances of these light elements with 
abundances that differ by nine orders of magnitude. 

There exist a variety of ways of modifying the standard BBN paradigm, in- 
cluding altering the spatial distribution of baryon number, out-of-equilibrium 
decays of massive particles, or new neutrino physics (for a summary, see 
Ref. [11]). 

We focus our attention on a minimal extension to BBN by a modification 
of the neutrino sector, and specifically to models which attempt to simulta- 
neously account for the indications of neutrino mixing and masses from the 
atmospheric neutrino results of Super-Kamiokande [12], the observations of 
the transformed solar neutrino flux [13,14], and the Liquid Scintillator Neu- 
trino Detector (LSND) signal [15]. A class of models that can accommodate all 
of these results introduce a fourth mass eigenstate [16]. As is well known, the 
fourth flavor state must be sterile, i.e., not participating in SU(2)l x U(1)y 
interactions, due to its being both light (m C 1 GeV) and not observed in the 
invisible width of the Z° boson [6]. In current manifestations of four neutrino 
mixing models, the sterile neutrino is not necessarily closely associated with 
a single mass eigenstate, since the atmospheric and solar observations each 
disfavor large sterile components. However, recent global analyses [16] of the 
available neutrino oscillation data, including short baseline limits [17], leave 
four-neutrino models viable. There are several ways that light sterile neutri- 
nos can be accommodated in neutrino mass models. For a review see, e.g., 
Ref. [18]. 

The standard contribution to the energy density by the three active neutrinos 
can be augmented by the complete or partial equilibration of the sterile mode. 
The increased energy density in units of the energy density in one neutrino and 
its antiparticle (p Pu = 7tt 2 T 4 /120) is then N v = 3 + AN u , where AA^ = p s /p Pu 
is the relative contribution of the sterile state. Here, we use AN U strictly as 
a parameterization of extra (or missing) relativistic energy density. Sterile 
neutrinos in the early universe can also give rise to lepton asymmetry genera- 
tion [19], which can alter or strongly suppress sterile neutrino thermalization, 
or, if the asymmetry is generated in the v e jv e sector, alter beta-equilibrium 
and thus light element abundance production, primarily in the production of 
4 He. As an alternate measure, the baryon-to-photon ratio rj is related to the 
cosmological baryon density Qb (as a fraction of the cosmological critical den- 
sity) as rj ~ 2.74 x I0~ 8 n b h 2 , where h is the present Hubble parameter in 
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units of 100 km s^ 1 Mpc -1 . 



Letting N v be a free parameter of BBN, its primary effect is altering the 
predicted 4 He abundance Y p . The remaining parameter, the baryon content 
77, is (over) constrained by D/H, 3 He and 7 Li. In one analysis, Lisi, Sarkar 
and Villante [20] used four permutations of primordial light element abun- 
dance determinations to derive limits roughly in the range 2 < N v < 4. In 
one combination of light element abundance determinations (their data set 
A), the 99.7% CL region allowed N v ~ 4.5. This value is widely cited as al- 
lowing for an additional neutrino (or relativistic degree of freedom) at BBN. 
Though this limit was correct, it relied on the possibility of a "high" primor- 
dial deuterium abundance. Since that work, deuterium has been observed or 
bounded to have a "low" value in six high-redshift quasar absorption systems 
(QAS) by three groups [9,21], and the "high" deuterium QAS observation [22] 
has not been verified in other systems and is disputed [23]. Using the low 
deuterium abundance and a small uncertainty in Y p , Buries et al. [24] found 
the limit N v < 3.20 with the prior N v > 3. The work by Cyburt, Fields & 
Olive [5] took the possibility of two values of the baryon density, that given 
by D/H+BBN and CMB, VL h h 2 ~ 0.02, and a value of Vt h h 2 ~ 0.01 preferred 
by one inferred primordial value of Y p [25] and the 7 Li abundance [26] (if un- 
depleted). Ref. [5] finds that these densities give 95% CL upper bounds of 
N v < 3.6 and N v < 3.9, respectively. 



Motivated by the six quasar absorption system measurements of a "low" D /H 
and the analysis of the observations of the CMB anisotropy experiments DASI, 
BOOMERanG, and MAXIMA, we adopt that the inferred value of the cos- 
mic baryon density from D/H plus standard BBN f^/i 2 (D/H) and the shape 
of the acoustic peaks in the CMB angular power spectrum fl b h 2 (CMB), are 
approaching the actual value of flbh 2 , within statistical and systematic uncer- 
tainties. The forthcoming analysis of the Microwave Anisotropy Probe (MAP) 
satellite's observations will potentially reduce the uncertainty in VL^h 2 to ap- 
proximately 10% [27]. In Section 2, we analyze in detail the constraints arising 
from accurate calculations of the primordial helium abundance, the inferred 
primordial helium abundance, using either f2 b /i 2 (D/H) or Of,/i 2 (CMB), and 
show that a thermalized fourth neutrino is highly disfavored by standard BBN. 
In Section 3, we analyze current four neutrino mixing schemes and their be- 
havior in the early universe and show that thermalization of a fourth neutrino 
state is unavoidable. In Section 4, we describe various means and methods of 
extending the standard model to evade this constraint. 
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2 The BBN Prediction of the Number of Neutrinos 



The consistency of BBN as a predictor of the light element abundances al- 
ready been explored in some detail [3,4,5,6]. We instead focus on the current 
uncertainties in the cosmic baryon density Q b h 2 and the observed primordial 
4 He abundance Y p , the light nuclide whose abundance is most sensitive to the 
energy content of the universe at the BBN epoch. 

The baryon content of the universe can be estimated in a variety of ways, of 
which the most precise measures currently are the deuterium abundance at 
high redshift and the shape of the acoustic peaks in the CMB [28] . Deuterium 
has been observed in high-redshift metal-poor neutral hydrogen systems which 
are seen as absorbers in the spectrum of back-lighting quasars. The deuterium 
in these extremely metal-poor systems is inferred to be close to the primor- 
dial value due to minimal stellar processing which produces metals and only 
destroys deuterium. Due to the extreme sensitivity of D/H to the baryon con- 
tent at BBN, the baryon density required to produce the observed deuterium 
abundance is rather precisely determined [4,9]: 

n b h 2 (D/E) = 0.020 ± 0.002 (95%CL), (1) 

for the standard energy density content N v = 3. The baryon density inferred 
from D/H has a small dependence on the relativistic energy density of the 
plasma, which is [29] 

ft b /i 2 (D/H,Ag = n b h 2 (D/R,N u = 3)(1 + 0.125 AN V ). (2) 



The baryon content also alters the amplitude of acoustic oscillations in the 
primordial plasma at CMB decoupling and the relative height of the first 
three acoustic peaks (for a review of the physics of the CMB, see Ref. [30]). 
The first three acoustic peaks in the angular power spectrum of the CMB 
have been detected in the analysis of CMB anisotropy measurements by the 
DASI [31], BOOMERanG [32] and MAXIMA [33] experiments. The results of 
these experiments' analyses find 

tt b h 2 (D) = 0.022±g;gg| (95%CL) 
Q b h 2 (B) = 0.022^;°o3 (95%CL) 
n b h 2 (M) = 0.033 ±0.013 (68%CL) (3) 

respectively. The BOOMERanG value above is that given by their Bayesian 
approach. For concreteness in our analysis, we quantify the uncertainty in the 
cosmic baryon density inferred from the CMB with the likelihood function 
given in Ref. [31] by DASI+DMR. 

If analyses of the CMB anisotropy measurements change and provide a value 
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for Qbh 2 that is higher than that inferred from standard BBN, then this could 
have been an indication for a model with large and disparate neutrino degen- 
eracy parameters known as degenerate BBN [34]. However, if the favored large 
mixing angle (LMA) neutrino mixing parameters of the solution to the solar 
neutrino problem are verified (e.g., by the KamLAND experiment [35]), then 
synchronized neutrino flavor transformation in the early universe stringently 
limits neutrino degeneracies [36] and degenerate BBN is no longer a rescue. 

In order to precisely predict the abundance of 4 He from standard BBN for 
varying baryon density and neutrino number, Lopez and Turner [37] included 
finite-temperature radiative, Coulomb and finite-nucleon-mass corrections to 
the weak rates; order-a quantum-electrodynamic correction to the plasma 
density, electron mass, and neutrino temperature; and incomplete neutrino 
decoupling. Ref. [37] provides a fitting formula for their results of the pre- 
dicted helium abundance as a function of 77, N u and neutron lifetime r. We 
employ this fitting formula for the predicted 4 He abundance, taking into ac- 
count the typographical correction of signs noted in Ref. [4]. After all of the 
above corrections are applied, the uncertainty in the predicted helium abun- 
dance is dominated by the neutron lifetime uncertainty, which is now known 
better than 0.1% [6]. Therefore, we can safely ignore the theoretical errors, as 
they are dwarfed by observational uncertainty, which we now address. 

The primordial helium abundance Y p has been estimated in observations of 
hydrogen and helium emission lines from regions of hot, ionized metal-poor gas 
in dwarf galaxies (Hll regions). By extrapolating the helium abundance and 
metallicity relationship for these regions to zero metallicity, Olive, Steigman 
and Skillman [25], and Fields and Olive [38] find 

F p (OSS-FO) = 0.238 ± 0.002 (stat.) ± 0.005 (sys.) , (4) 

while Izotov and Thuan [39] find 

Y p (IT) = 0.244 ± 0.002 (stat.) . (5) 

Uncertainties regarding the ionization structure and temperature uniformity of 
the Hll regions as well as underlying stellar absorption are sources of significant 
systematic error. Refs. [8,25] estimate systematic effects in the primordial 
helium abundance are 2%. Ref. [7] finds that systematic effects can lead to 
2-4% uncertainties that tend to overestimate Y p . In an attempt to avoid bias, 
in this work we adopt the central value of 

Y p = 0.241 ± 0.002 (stat.) ± a sys , (6) 

and characterize the systematic uncertainty as the disparity between compet- 
ing claims 

a sys = |y p (OSS-FO)-y p (IT)|, (7) 
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or approximately 3%. The shape of systematic uncertainty in likelihood space 
is certainly not well defined, therefore we make the simplifying ansatz of a 
Gaussian distribution, as done, e.g., in Ref. [5], and combine the statistical 
and systematic errors in quadrature. 

We produce probability distribution functions (p.d.f.'s) for N v versus Q b h 2 , 
using Gaussian distributions for Y p [Eq. (6)] and Q b h 2 (D/B., N u ) [Eq. (2)], and 
the likelihood function given in Ref. [31] for Q b h 2 (T)ASl). We find, using either 
the information from deuterium or the CMB on VL b h 2 : 



JV,(D/H) = 2.60l°J° (8) 
JV„(DASI) = 2.461*;°?, (9) 

at 95% CL, which is consistent with the standard BBN prediction. We show 
the shapes of the likelihood contours in Fig. 1 (a). To illustrate the difference 
between adopting the IT or OSS-FO helium values, we plot the 99% likelihood 
contours for the choices Yp(OSS-FO), Y p (lT) (using the systematic uncertainty 
of ±0.005) and our choice (6). As seen in Fig. 1 (b) the range of uncertainty 
does not depend on the choice of the central value but the size of systematic 
effects. 

For a sterile neutrino to be thermalized with the bath of the early universe, 
the active neutrinos must be thermalized initially. This constitutes prior in- 
formation that may loosen the constraints shown in Fig. 1. Prior information 
can be included in a Bayesian approach [6], integrating the p.d.f. only in the 
physically allowed region, N u > 3, which we have done using a Monte Carlo 
integration. We show the confidence level intervals for this case in Fig. 2. As 
seen there, a fully populated fourth neutrino is excluded at approximately the 
Yo CL. 



Measurements of CMB anisotropics by the MAP satellite may measure VL b h 2 
to 10%, giving a precise value independent of BBN. If consistent with Q b h 2 
inferred from D/H and BBN, then Q b h 2 becomes a "nuisance parameter" that 
can be marginalized and the confidence level for N u is simply the integral of 
the p.d.f., 

CL(N„)= p(N' u \Y p )dK. (10) 

J3 



3 Four Neutrino Models in the Early Universe 



There is now convincing evidence for neutrino flavor states to be composed of 
large amplitudes of more than one mass state from two experiments: Super- 
Kamiokande [12] and the Sudbury Neutrino Observatory (SNO) [14]. There 
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Fig. 1. Shown are the contours of 68%, 95% and 99% CL (inner to outer con- 
tours) when using the baryon density inferred from the deuterium abundance 
(D/H) at high-redshift (solid lines) or that from the DASI+DMR analysis of CMB 
anisotropies (dashed lines), in frame (a). In both cases, we use our our adopted 
determination of the primordial 4 He abundance. In frame (b), we show the 99% CL 
contours using the Izotov & Thuan (IT), Olive, Steigman & Skillman (OSS) and 
Fields & Olive (FO) and our adopted value and uncertainty of the primordial 4 He 
abundance along with the D/H determination of S7f,/i 2 . The uncertainty in N u does 
not depend on the choice between IT and OSS-FO as much as the size of systematic 
uncertainty. See text for details. 



is also an indication of a neutrino oscillation signal at short baselines from 
the Liquid Scintillator Neutrino Detector (LSND) experiment [15]. To accom- 
modate all three of these results, a four neutrino model must be invoked (or 
CPT is violated; see below). The mass and flavor state bases are related by a 
unitary transformation 

4 

v a = ^2U ai Vi, (11) 

i 

where a = e, /x, r, s denotes the flavor state and i is the mass state. The matrix 
U a i generally has 6 rotation (mixing) angles and 3 CP violating phases. The 
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Fig. 2. Shown are the contours of 68%, 95% and 99% CL (inner to outer contours) 
given the prior condition that the active neutrinos are thermalized, our adopted 
value of primordial 4 He abundance, and f2&/i 2 determined by D/H (solid lines) or 
the CMB observations by DASI+DMR (dashed lines). See text for details. 

transformation probability has the form: 



P(v n 



Ja/3 



^^(U^U^U*^) sin 2 (d 



i>3 



+ 2 £ Im (u^UpiU^) sin Um\ 



i>3 



13 IE) ' 



'121 



which provides rich neutrino oscillation phenomenology and experimental pos- 
sibilities [40]. Here Smfj = m 2 — m 2 , L is the distance from where the flavor 
u a was created, and E is its energy. For reviews of neutrino phenomenology 
see, e.g., Refs. [41]. 

For neutrino physics in the early universe, we are interested in the magnitude 
of mixing amplitudes of active neutrinos converting to sterile states. First 
derived by Langacker [42], and Dolgov and Barbieri [43], it has been known 
for some time that a sterile neutrino coupling with a single active neutrino via 
the unitary neutrino mass matrix must not have [44,45] 



5m as sin 2#bbn < < g ^ 1Q _ 6 



for a = e 
for a = (jl,t 



(13) 



in order to not fully thermalize the sterile neutrino prior to BBN via non- 
resonant collisional processes. Such constraints (13) certainly do not directly 
apply to multiple neutrino mixing schemes including a sterile which nature 
may have given us. Multiple mixing angles and the phenomenon of lepton 
number generation via neutrino mixing complicate the BBN bound. 



We adopt a rotation angle ordering for U so that [y e , v^, v T , v s ) = z/ 2 , z/ 3 , i/ 4 ) 
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when all mixing angles are set to zero. Four neutrino mixing scheme constraints 
have been examined recently in detail by Di Bari [46], which we summarize 
and expand on here in view of the recent global analyses of four-neutrino 
models by Maltoni, Schwetz and Valle [16] . The effective oscillation amplitude 
of active-sterile neutrino mixing between flavor a and the sterile can be written 
as 

A a . s = A\U a4 \ 2 \U s4 \ 2 ~ sin 2 2#bbn • (14) 

Therefore, BBN constraints for active-sterile mixing through a pair of neutrino 
mass eigenstates z/4, z/j are 

x 2 A 2 <r J 5 x 10 ~ 6 > for a = e 

Sm^A <l , (15) 

13 x 10 , for a = /i, r 

for nonresonant sterile production, where Sm^ > 0. This constraint (15) ap- 
plies primarily to the active-sterile mixing that leads the oscillation, i.e., that 
which has the shortest oscillation length or largest 5m 2 -. Small vacuum mix- 
ing amplitudes in the leading oscillation mode may avoid thermalization, and 
so secondary oscillation modes with larger mixing amplitudes may thermalize 
the sterile. 



3.1 3+1 



Models referred to as (3+1) may satisfy all experimental indications of neu- 
trino oscillations with a triplet of mass eigenstates that provide the atmo- 
spheric and solar mass-scales, and a sterile-dominated mass eigenstate with 
a large mass-scale splitting with the triplet providing the LSND result via 
indirect — > v e mixing through the sterile. In order to satisfy the mixing am- 
plitude that would provide the LSND signal, the oscillation amplitude must 
be [16] 



Afj,- e = A\U e4 \ 2 \U^\ 2 

> 3 x 1(T 4 (99% CL) , (16) 

at a 5m1 SND — 2 eV 2 from Fig. 8 of Ref. [16]. Indirect mixing of this form has 
two mixing amplitudes that may thermalize the sterile. In this case, the mixing 
amplitudes A^ s and A e . s may participate in the thermalization. Consider the 
slightly less-constrained [cf. (15)] A^ ;s = A\U^\ 2 \U si \ 2 . The mixing matrix 
element 

\U sA \ 2 > 0.54 (99% CL) , (17) 

is bounded from below from constraints on the fraction of sterile neutrinos 
participating in atmospheric oscillations (see Fig. 3 of Ref. [16]). 
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Evading constraints from BBN on A^ s and A e - S would require minimizing both 
\U e ±\ and \U^\ while satisfying (16). This gives |t/ e4 | 2 = \U^\ 2 ~ 10~ 2 . Com- 
bined with the limit (17), the BBN constrained combination has the minimum 
value 

Smt SNB A% s > 7 x 10- 4 (18) 

which exceeds the limits (15) by at least two orders of magnitude and in- 
variably thermalizes the sterile. This constraint comes from the conservative 
case where Sm^ > 0. The inverted case Sm^ < is resonant and more strin- 
gently constrained. Therefore, (3+1) models are strongly disfavored by stan- 
dard BBN. 



3.2 2+2 



Four neutrino models may also accommodate all indications for neutrino mix- 
ing with mass eigenstates in a pair of doublets that provide the solar and 
atmospheric mass scales and a large mass gap between the doublets providing 
the LSND mass scale. The global analysis by Maltoni et al. [16] finds that such 
(2+2) models are consistent within 99% CL either with the sterile neutrino 
completely participating in the atmospheric or solar solutions. Therefore, it is 
possible to choose a very small amplitude mixing between the doublet not par- 
ticipating in sterile oscillations and the sterile state such that the large LSND 
mass splitting does not populate the sterile neutrino. On the other hand, uni- 
tarity constrains the sterile state to be present among some linear combination 
of mass eigenstates. Whether (2+2) scenarios are compatible with the exclu- 
sion of large sterile components in both the atmospheric and solar neutrino 
observations is controversial [47]. 

Whether the sterile flavor is in the atmospheric doublet or solar doublet is 
not of concern for the early universe. In either case, the sterile flavor will be 
thermalized. For our notation, we employ the mass scheme used in Ref. [16], 
where the (2+2) model has the solar scale between v\ and z/4 and the atmo- 
spheric scale between z/ 2 and z/ 3 , with the hierarchy only being determined by 
the condition for resonance in the sun, m x < m 4 . One can avoid the effects of 
thermalization of the largest mass scale by setting the inter-doublet mixings 
to zero, U s i = U s ^ = 0, but by unitarity \U S 2\ 2 + \U S i\ 2 = 1, whereby the 
sterile neutrino participates in large part in the solar scale, atmospheric scale 
or both. Having a complete sterile solution for either scale has already been 
known to thermalize the sterile neutrino [44,45]. One could consider demo- 
cratically separating the sterile into both the atmospheric and solar scales to 
minimize its presence in both, so that \U S 2\ 2 = \U S i\ 2 = 1/2. However, the os- 
cillation amplitudes A^g = 4:\U /1 2\ 2 \U S 2\ 2 and A e - S = 4|f/ e 4| 2 |f/ S 4| 2 still grossly 
exceed the BBN bounds (15) since the magnitudes |?7 e 4| 2 = \U e \\ 2 tan 2 6*lma 
and I 1 must be large to accommodate the large to maximal mixing angle 
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solutions of the solar and atmospheric neutrino problems. Therefore, (2+2) 
models are also strongly disfavored by standard BBN. 

3.3 Self-Suppression 

The possibility that a four neutrino mass scheme could be arranged in such 
a way as to evade sterile-thermalization constraints were considered by Bell, 
Foot & Volkas [48] and Shi, Fuller & Abazajian [49]. One could potentially 
either self-generate a lepton number and suppress the large-mixing-amplitude 
thermalization or offset the effects of sterile thermalization by altering the elec- 
tron neutrino-antineutrino asymmetry through alteration of beta-equilibrium, 

n + u e <-> p + e~ 

n + e + ^p + u e . (19) 

In the models considered in Refs. [48,49], the direct thermalization bounds (15) 
were avoided by placing the sterile neutrino in the small-mixing-angle solution 
to the solar neutrino problem, a region of parameter space still viable at the 
time and outside of the constraint region (15). 

In addition, Refs. [48,49] explored methods of generating asymmetries be- 
tween electron neutrinos and antineutrinos by resonant lepton number gener- 
ation [19]. The resonance condition in the early universe requires m 4 < rrii, 
where rrii is a mass eigenstate (more) closely associated with an active flavor. 
A positive electron neutrino number will suppress the 4 He abundance by shift- 
ing the rates (19), which would be necessary if standard BBN is inconsistent 
by having too high of a predicted 4 He abundance for a given Vt b h 2 . As shown 
above, standard BBN remains consistent within observational uncertainty. The 
sign of resonantly generated electron neutrino/antineutrino asymmetry can be 
chaotic [50], or at least not well determined [51], having an significant chance 
(50% if chaotically random) of being negative and actually increasing Y p by 
altering beta-equilibrium in the opposite direction. If the sign of the asymme- 
try is randomly chaotic, then causally disconnected regions will have different 
sign asymmetries, which leads to an enhancement of the transformation of 
active neutrinos into sterile neutrinos at the boundaries of regions of different 
sign [52] and potentially placing more stringent constraints on four neutrino 
mass schemes [53]. 

There is considerable evidence now that the solar solution lies in the LMA re- 
gion of parameter space [14]. Therefore, as discussed in the previous sections, 
thermalization of the sterile is unavoidable in either the (3+1) or (2+2) sce- 
narios. And, importantly, it was shown by Di Bari [46] that thermalization of 
the sterile in these four neutrino models suppresses lepton number generation, 
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and electron neutrino/antineutrino asymmetries are not effective in avoiding 
the BBN bounds (15). 



4 Constraint Evasion and New Physics 

The simplifying and appealing principle of Occam's razor has proven to be a 
powerful tool as a predictor in science, yet nature does not always take the 
most simple form. The minimal model for four neutrino mixing or the standard 
BBN described above may certainly not be the entire framework of the early 
universe or particle physics. Importantly, if all experimental indications for 
neutrino oscillations remain, viz., if the MiniBooNE detector [54] verifies the 
LSND signal, K2K [55] and MINOS [56] verify the atmospheric oscillation 
solution and KamLAND detects the LMA signal [35], then new physics must 
be at play beyond standard three-neutrino mixing and standard BBN. There 
exist a number of ways of accommodating such a scenario, several of which are 
described below. The aesthetic value of these scenarios are left to the judgment 
of the reader. 

Pre-existing lepton asymmetry — A lepton number in the active neutrino 
flavors will suppress sterile neutrino population by magnifying the associated 
lepton potential and dwarfing the vacuum mixing amplitude [57] . This lepton 
number would have to be produced by an unspecified mechanism earlier than 
the population of the sterile neutrino would take place. 

A fifth mass eigenstate — Appropriate insertion of a mass eigenstate with 
a major sterile component with m 5 < rrii (i = 1...4) in degenerate neutrino 
mass models may resonantly generate lepton number sufficiently prior to sterile 
thermalization as to suppress it. This possibility was explored in Ref. [46]. 

Majoron fields — One mechanism for generating neutrino mass involves a 
massless Nambu-Goldstone boson (a majoron) from models where either the 
total or partial lepton number is spontaneously broken [58]. In such models, 
a coherent majoron field creates potentials for the neutrinos proportional to 
the gradient of the field and suppresses sterile thermalization in a similar way 
as a pre-existing lepton asymmetry [59]. Interestingly, this mechanism arises 
from the neutrino mass model itself. 

A low reheating temperature universe — There is no direct evidence that the 
neutrino background is thermalized. As explored in Refs. [60], the highest 
temperature of the universe could have been only 0.7 MeV. The neutrino 
background may never have been thermalized, but the observed light element 
abundances could still be created. In this case, sterile neutrinos may modify the 
nucleosynthesis processes by partial population but are not directly excluded. 
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Baryon-antibaryon inhomogeneity — Detailed calculations of diffusion and 
nucleosynthesis in universes containing baryon number asymmetries [61] have 
found that small-scale antibaryon domains are not excluded by BBN and the 
observed light element abundances [62] , and may lift constraints on relativistic 
energy density present at BBN to N v < 7, even with total baryon densities 
consistent with CMB observations [63]. 

Extended quintessence — Non-minimally coupled quintessence models (where 
the quintessence field is not only coupled to gravity) that provide a negative- 
pressure vacuum energy density to explain the acceleration phase that the 
universe may be entering can alter BBN [64]. In certain cases of such "extended 
quintessence" scenarios, the quintessence field may behave to decrease the 
expansion rate during the freeze-out of beta equilibrium (19), and therefore 
decrease the predicted helium abundance. This reduction of the expansion rate 
could offset the increase in the expansion rate due to the presence of an extra 
neutrino degree of freedom and allow for four-neutrino models. 

CPT violating neutrinos — There exists a radical proposal that fits all in- 
dications for neutrino oscillation and invokes CPT violation in the neutrino 
sector [65]. The success of this model lies in the fact that LSND's indication 
for neutrino oscillation lies primarily in the antineutrino P M — > v e channel [15], 
is motivated by braneworld scenarios with extra dimensions, and gives dra- 
matic predictions for the MiniBooNE [54] and KamLAND [35] experiments. 
This model has no effect on standard BBN since in the standard case the neu- 
trinos and antineutrinos are equally thermally populated and therefore CPT 
violating neutrino oscillations do not disturb the detailed balance of thermal 
equilibrium, leading to no direct conflicts between light element abundances 
and standard BBN. 



5 Discussion and Conclusions 

In minimal models of big bang nucleosynthesis with no new physics, we have 
shown that four neutrino models explaining current indications for neutrino 
oscillations are disfavored at the 99% CL. This conclusion depends on system- 
atic effects not being larger than that expected (~ 3%) in determining the 
4 He abundance in ionized Hn regions and that the baryon density inferred by 
the D/H abundance in six high-redshift quasar absorption systems and the 
anisotropies in the cosmic microwave background are approaching the true 
cosmic value of VL b h 2 . 

The MAP satellite will verify or disprove the value of VL b h 2 inferred by the 
above methods to high precision in the near future [27]. If consistency re- 
mains in cosmological determinations of the baryon density, then Q b h 2 would 
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then become a "nuisance" parameter in determining the cosmological energy 
density at standard BBN. The primordial helium abundance is currently the 
best probe of the energy density of the universe present at BBN, yet there 
remain no concrete proposals in the literature for the reducing systematic 
uncertainties present in determining the primordial 4 He abundance, which is 
the dominant uncertainty in constraining the energy density present in the 
universe at the age of one second via standard BBN. 

In addition, we have summarized several scenarios that evade the standard 
BBN model constraints presented here. Remarkably, if all experimental indi- 
cations for neutrino oscillations are confirmed, new physics must be present 
not only in the particle content of the neutrino sector but also in the early 
universe. 
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